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Continuous measurements play a pivotal role in the study of dynamical open quantum systems. 
‘Dyne’ detections are among the most widespread and efficient measurement schemes, and give rise 
to quantum diffusion of the conditioned state. In this work we study under what conditions the 
detector dependency of the conditional state of a quantum system subject to diffusive monitoring 
can be demonstrated experimentally, in the sense of ruling our any detector-independent pure-state 
dynamical model for the system. We consider an arbitrary number L of environments to which 
the system is coupled, and an arbitrary number K of different types of dyne detections. We prove 
that non-trivial necessary conditions for such a demonstration can be determined efficiently by semi- 
definite programming. To determine sufficient conditions, different physical environmental couplings 
and Hamiltonians for a qubit, and different sets of diffusive monitorings are scrutinized. We compare 
the threshold efficiencies that are sufficient in the various cases, as well as cases previously considered 
in the literature, to suggest the most feasible experimental options. 
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I. INTRODUCTION 

One of the most challenging tasks in the analysis of an 
open quantum system is to gain as much information as 
possible about the system state. The conventional ap¬ 
proach of tracing over the environment, with which the 
system is coupled to, would not yield to the maximal 
knowledge. However, continuous measurements (also re¬ 
ferred to as monitorings ) of the environment make it pos¬ 
sible to gain further information about the system by ex¬ 
ploiting quantum correlations between the two [I]. The 
former approach results in a master equation (ME) for 
the state matrix p of the system whereas the latter gives a 
stochastic master equation (SME) describing evolution of 
a conditional state matrix g given a readout. The latter 
is known as an unraveling of the ME [ 2 ] , because the ME 
is recovered by averaging over all possible measurement 
records. Under ideal conditions, for every single mea¬ 
surement record there is a path of a ray in the system’s 
Hilbert space, called a quantum trajectory. 

A SME may have more than one physical realization 
(that is, the set of different ways of monitoring the envi¬ 
ronment is a surjective mapping onto the set of all pos¬ 
sible SMEs). In the context of quantum optics, differ¬ 
ent monitorings of the environment can be achieved by 
adding a local oscillator (LO) to the system output field 
prior to detection. Depending on the strength, frequency, 
and phase of the LO in comparison with those of the emit¬ 
ted radiation the dynamics of the system are described 
by distinct SMEs. These fall into two categories: jump 
unravelings and diffusive unravelings |pQ. For quantum 
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jumps the LO either is a weak field or is set at zero (di¬ 
rect detection), while for quantum diffusion it is a strong 
field, which can be resonant with the output field (ho¬ 
modyne detection) or detuned from it (heterodyne de¬ 
tection) |TJ[3j. 

In this paper we scrutinize diffusive unravelings of MEs 
in the context of experimentally demonstrating the de¬ 
tector dependency of quantum dynamics [1]. Historically, 
before the emergence of quantum trajectory theory 0 ca¬ 
lm, it was widely believed that a system’s quantum dy¬ 
namics was objective, and thus independent of any mea¬ 
surement. This notion was proposed in early works of 
Einstein on emission and absorption of light quanta by 
an atom [9] such that the system (atom) state under¬ 
goes evolution among stationary states of Bohr [TO]. But 
since the early 90’s, modern theory of quantum dynamics 
has rendered this idea obsolete, implying that dynamical 
quantum events are subject to a remote detector. 

It remains a challenge, however, to verify experimen¬ 
tally the testimony of this theory. In 2012, a proposal 
was made for an experiment to disprove all objective 
pure-state dynamical models (OPDM) for a particular 
quantum system [4] . If an experimentalist had per¬ 
fect detection (and collection) efficiency, r) = 1 , the con¬ 
ditioned pure state would be pure, and it would sim¬ 
ply be a matter of verifying that two different detection 
schemes lead to different types of pure states. That is, 
the problem would be a continuous-in-time version of the 
Einstein-Podolsky-Rosen (EPR) phenomenon [XT] . How¬ 
ever, in the real experimental setups 77 < 1 , and in or¬ 
der to rule out all OPDMs one would require to invoke 
the more general form of the EPR argument known as 
steering [T2] (the term introduced by Schrodinger for this 
phenomenon [13]), or EPR-steering [14]. In this formal¬ 
ism one experimenter (Alice) proves that the objective 
quantum state of the other experimenter (Bob) cannot 
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be explained by any local hidden state (LHS) models. 
This is fulfilled by expressing the correlation of measure¬ 
ments of the two parties in terms of an inequality, the vi¬ 
olation of which can be interpreted as ruling out all LHS 
models for Bob’s system. Generalizing this to continuous 
monitoring [4], inequalities can be derived of which the 
violation would rule out all OPDMs, and consequently 
confirm the detector-dependency of stochastic quantum 
dynamical evolutions. 

The physical system studied in Ref. [4j was a reso¬ 
nantly driven two-level atom with two measurement ar¬ 
rangements to capture the fluorescence. It turns out 
that using two homodyne detection, and assuming both 
have equal efficiencies, the critical efficiency to disprove 
OPDMs is rj c ss 73%, which is quite high. Further, it was 
proved H] that p > 50% is a prerequisite for being able to 
demonstrate EPR-steerability using any dyne detection 
schemes. Very recently, two of us have generalized this 
necessary condition [TjDj, showing that for any Marko¬ 
vian open quantum system with an arbitrary number L 
of decoherence channels to the environment and with an 
arbitrary number K of diffusive unravelings, unless at 
least one member of the set {rj^} is greater than 50%, 
the diffusive unravelings cannot be used to demonstrate 
detector dependency of the conditional state matrix g. 

The present work furthers the investigation of how well 
the detector dependency of dynamical quantum events 
can be proven experimentally using diffusive unravelings, 
in terms of the efficiency required. Our first aim is to 
derive even more general necessary conditions, and our 
second aim is to derive sufficient conditions in a variety 
of scenarios with a simple system (a qubit). 

We organize the paper as follows. First, Sec. [n] re¬ 
views briefly formalism of general diffusive unravelings. 
In Sec. m it is then shown that the most general proof 
of a no-go for inefficient diffusion can be cast as a spe¬ 
cial case of a semidehnite programming (SDP) known as 
a feasibility problem (FP). We present a simple graphi¬ 
cal representation of this general necessary condition for 
the case where the system interacts with just one envi¬ 
ronment. Following this, In Sec. |IV| different MEs and 
measurement strategies are considered. First we ana¬ 
lyze a system with three irreversible channels, and three 
different diffusive unravelings to try to prove detector- 
dependence. Next, systems with fewer decoherence chan¬ 
nels and just two different diffusive unravelings are stud¬ 
ied. Finally, in Sec. [V] we conclude by summarizing the 
results of this study along with the outcomes of previous 
publications in Table [T| 

II. GENERAL DIFFUSIVE UNRAVELING 
WITH IMPERFECT DETECTION 


Here H is the system Hamiltonian (in units where H = 1), 
c = (ci,.. ., c l ) t is a vector of arbitrary operators (also 
know as Lindblad operators PI]), and V[c] = Ya = 1 v [^i] 
where the decoherence superoperator V[d\p = opo' — 
1/2{<5^6, p}, with curly brackets representing the anti¬ 
commutator. 

The most general diffusive unraveling of the ME, 
Eq. {!]), by allowing for inefficient detection reads [I] 

dg = Cgdt + U[dV^c}e, (2) 

where T-L[d\p = dp + pd ^ — T r[op + pd^]p is a nonlinear 
superoperator, and dV = (dVi ,..., ciVi) T is a vector 
of infinitesimal c-number Wiener increments. These are 
Gaussian random variables satisfying E[dV] = 0, where 
E[o] denotes the ensemble average of the random variable 
o with respect to its probability distribution function. 
Physically, they should actually be regarded as noise in 
the output complex photocurrent 

J dt = (0 c + T c*) dt + dV, (3) 

with the following correlation relations 

dVdV f = Odt, dVdV T = T dt, (4) 

where = (o T )^, and (6) = Tr[dp] is the quantum 
mechanical expectation value of an operator 6 with re¬ 
spect to the conditioned quantum state g. Here 0 = 
diag(? 7 i,..., pl) is a real diagonal matrix allowing for 
imperfect detection in the formalism so that channel l 
is monitored with efficiency 0 < rji < 1. The rest of pa¬ 
rameters that characterize diffusive unravelings are being 
encoded in a complex symmetric matrix T = Y T . It is 
convenient to encapsulate all of the correlation properties 
in the so-called unraveling matrix na 


C/(0, T) 


1 / 0 + Re [T] Im [T] \ 

2 ^ Im [T] 0 — Re [T] ) ' 


( 5 ) 


with the only constraint that it must be positive semi- 
definite (PSD), i.e. that 3 Z £ £ 2 Lx 2 L g^p that 
[7(0, T) = Z T Z. Of course, one could work with other 
matrix representations of diffusive measurements other 
than U 118| , but for the sake of consistency with Ref. [T5] 
we rather keep this notation, even though we shall also 
give a technical reason why it is not an appropriate choice 
(see Sec. Ill B). Thus an ideal monitoring of the envi¬ 
ronment giving rise to a pure-state quantum diffusion is 
represented by U(I, T), where I £ 1R L is an identity ma¬ 
trix, and correspondingly Eq. ([2]) can be replaced with a 
stochastic Schrdinger equation. 

An interesting class of unravelings is obtained when 
the Lindblad operators are Hermitian c = and the 
complex Wiener increments satisfy 


The state p of an open quantum system in the Markov ( ]y l = e^^/rjdWi, l = 1, • • • , L, 

approximation, defined by tracing over the environment, 
has its evolution governed by a ME where 

p = -i [H, p] +V[c] p = Cp. 


( 6 ) 


( 1 ) 


dWj dWk = Sjk dt. 


( 7 ) 
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Therefore, depending on the monitoring scheme the 4h 
measurement outcome would be either a current contain¬ 
ing maximal information about the observable c/> 

Jf UIient dt = 2r ] (c l )dt + ^TjdWi (0 = 0), (8) 

or just a pure-noise, 

Jf oise dt = iy/rjdWi (0 = tt/2). (9) 

In the first case, the associated conditional evolution 
tends to localize the system to a cj eigenstate: 

L 

dg = £gdt + y/rj ^ [(cj - (ci))g + R.c.]dWi. (10) 
1=1 

In the latter case, the stochastic evolution corresponds 
to a noisy Hamiltonian, since no information about the 
system is being obtained: 


dg = Cg dt + iyfr\ E [eA]dW t . (li) 


i=i 


We employ these kind of unravelings in Sec. IV for some 
specific scenarios. 

Now having this formalism in our arsenal we can inves¬ 
tigate proving detector dependency of diffusive stochastic 
evolution. 


III. NECESSARY CONDITIONS FOR 
DEMONSTRATING DETECTOR DEPENDENCE 

A. Coarse graining 

In addition to the mathematical tools that have been 
described thus far we require another essential concept, 
that is, coarse graining of diffusive unravelings. To this 
end, consider two unravelings U = 1/(0, T) and Uq = 
I/(0o, T 0 ) with associated vectors of Wiener processes 
dV and dV o, respectively. If it is possible to express the 
latter as the algebraic sum of the former and another 
complex vector Wiener increment dV obeying the above 
restrictions, i.e. dVo = dV+dV , then the first unraveling 
U can be realized experimentally by implementing the 
second one, Uq, and keeping just the relevant information 
in record, discarding the information in dV . This will be 
so if the tacit unraveling matrix for dV, U = Uq — [/, 
is also PSD. Under these conditions we call U a coarse 
graining of Uq, and Uq a fine graining of U (and of U). 


B. Arbitrary number L of environments 

Consider an experimenter who is able to implement a 
set it of some number K > 1 of distinct unravelings 


As pointed out in the introductory section, it is necessary 
to perform EPR-steering in order to prove that condi¬ 
tional stochastic evolution of the system is determined by 
a distant detection apparatus. This means that if there 
exists a single unraveling Uq so that V k, Uk = Uq — Uk € 
ip 2 Lx 2 L pgp> then every member of the set il is a coarse 
graining of Uq. If such a Uq exists, then 0 O can always 
be chosen to equal I (as this just makes every Uk more 
positive), so that Uq is a purity-preserving unraveling. In 
other words, there exists a pure-state dynamical model 
which could underly all of the different stochastic evolu¬ 
tions induced by the K different measurement schemes, 
namely that corresponding to Uq. Consequently, that set 
il cannot possibly be used to rule out all OPDMs. 

Thus we have a necessary condition for demonstrat¬ 
ing detector dependence: the non-existence of a common 
fine-graining Uq for every member of il. But is there a 
way to establish the existence, or non-existence, of such 
&Uq given it? 

In Ref. [15] a partial answer was given, in the form 
of a no-go theorem for inefficient diffusive unravelings: 
if for all decoherence channels l and for all monitoring 
schemes k the detection efficiencies satisfy rtf < 0.5, then 
there exists an unraveling Uq from which the set il can 
be obtained by coarse graining. This was proven by con¬ 
sidering the choice Uq = U(I,0). This unraveling cor¬ 
responds to “quantum state diffusion” (QSD) as intro¬ 
duced in Refs. [T9] [20j (as an objective pure state dy¬ 
namical model) and could be realized by unit-efficiency 
heterodyne detection [5j. 

In this paper we prove a much more general no-go the¬ 
orem. Rather than considering a particular Uq, and see¬ 
ing what that implies, we consider any given set 11 and 
define a procedure to determine whether any Uq exists. 
The procedure is, in general, numerical rather than an¬ 
alytical, but it can be performed efficiently because of 
the power of semi-definite programming [22] - To reiter¬ 
ate the problem definition, given a set of diffusive un¬ 
ravelings if = {t/(0fc, Tfc)}]^, we want to determine 
whether there exists a To such that [/(/, To) > 0 and 
Vfc,J7(J-0 fc , T 0 —T fc )>0. 


Solution via Semi-Definite Programming 

With this formulation in hand, the problem can be 
formalised as a particular class of semidefinite program¬ 
ming (SDP) which is referred to as a feasibility problem 
[21], the description of which comes below. The standard 
definition of a (dual) semidefinite program is |22j 

minimize b T x (13) 

subject to F(x) > 0, (14) 

where 


F(x)=F° + Y,XjF i - 

/=i 


11 = {14 = J7(0 fc , T fc ) : U k > 0, k = (1,..., K)}. (12) 


(15) 
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Here b £ R n , and Fq, ■ ■ ■ ,F n £ IR mxm are the known 
data, and the minimization is over x £ IR”. In terms 
of these variables, the constraint (14) is a linear matrix 
inequality (LMI). In the special case in which b = 0, 
the optimization problem reduces to the search for some 
vector x satisfying the LMI, which should obviously be an 
easier computational job in practice. This is an instance 
of a feasibility problem and is exactly what we need. 

To apply the above formalism, we need a one-to-one 
linear mapping between the symmetric complex L x L 
matrix Tq and a real vector x £ (R i ( i+1 ) [23] 



To = <?(x), 


( 16 ) 


Defining 4|l(x) = [/(/,^?(x)), which is still a linear func¬ 
tion of x, the required constraints can be written as the 
LMI (14) by choosing F(x) to be the following block- 
diagonal matrix £ S 2i ( if + 1 ) : 


F(x) = 


/ *(x) 

0 


0 ^ 

0 

♦ (x) - Ui 


0 




0 

\ 0 

0 


♦ (x) - U k j 


■ (17) 


This is because positive semidefiniteness of Eq. en¬ 
sures that all submatrices must be PSD. Also the linear¬ 
ity of <|k(x) and consequently that of E(x) implies that 
the latter can be written in the form of Eq. (15), with 
F° ,F n £ S 2L ^ +1 \ where n = L{L + 1). Therefore 
the task is to find the following feasible set 


FIG. 1. Feasibility set and space of unravelings for an open 
quantum system with one environmental decoherence channel 
L = 1. The light gray cones represent all physical unravelings 
that satisfy Eq. (201, and the dark gray cones show unrav- 

for (a) a 


elings that meet the constraints given in Eq. ( |21b| ) 
generic unraveling Uo = U(jio,vo), and (b) the QSD unravel¬ 
ing Uo = 17(1, 0). The intersection of these two conical spaces 
encloses a volume so that all points inside of it can be a coarse 
graining of the Uo located at the vertex. 


C. Single environment L = 1 

There are numerous physical problems that are mod¬ 
elled using the ME given in Eq. ([I]) with only one deco¬ 
herence channel, such as a two-level system with dipole¬ 
coupling to a coherent electromagnetic field [26]. In this 
case the vector of arbitrary operators c in Eq. |2| is re¬ 
placed with a single Lindblad operator c and hence the 
diffusive unraveling is described by 


U(V, T) = 


1 


r] + Re [w] 
Irn M 


Im [u] 

77 — Re [u] 


(19) 


S'= {x G R i(i+1) : F(x) > 0}. (18) 

Each member x of this set is called a feasible point , 
and defines a fine graining 4|k(x) of every member of the 
sets if as coarse grainings of it. This FP has 2 L[K + 1) 
constraints with L(L+ 1) unknown real parameters. This 
suggests that other representations of diffusive unravel¬ 
ings which deal with a larger number of unknowns, for 
example the one introduced in Ref. m , are not com¬ 
putationally efficient for the particular kind of problems 
we are considering here. It is known that a SDP can 
be solved in (approximately) polynomial time with re¬ 
spect ton = 2 L(K + 1), and with errors as small as one 
wishes [24] . Thus, finding the unknowns would not be 
costly from computational viewpoints [25] . in particular 
because no optimization is performed. If it turns out that 
the set is null then all OPDMs will be ruled out. 

This is the no go theorem. It is easiest to see how 
it works by considering the simple case L = 1; that is, 
when the system is allowed to couple with just one envi¬ 
ronment. 


where 77 G IR and v £ L are scalars corresponding to ma¬ 
trices 0 and T in Eq. ([ 5 ]), respectively. It is straightfor¬ 
ward to check out the positive semidefiniteness constraint 
for Eq. (19) by inspecting its spectrum. This constraint 
translate to 


M < v- 


( 20 ) 


Also, the same procedure applies to the K +1 constraints 


of the feasibility set, Eq. (18), so that they all can be 


arranged in the following form 


M < Vo, 

\/k \v 0 - v k \ < rjo ~ Vk- 


( 21 a) 

( 21 b) 


Since each diffusive unraveling can be uniquely deter¬ 
mined by (L 2 + 2L)|y = i = 3 variables, it is then possible 
to visualize a feasible set in three dimensional space such 
that the coordinates are labeled by Re[y],Im[u], and 77 . 
The set of all physically meaningful unravelings that sat¬ 
isfy Eq. (20) is shown in Fig. la as the light gray cone. 


Therefore the set it given by Eq. (12) that is supposed 


to be implemented by an experimenter is represented by 
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a set of points ffc)}fLi i n this cone. Moreover, it is 
evident that unravelings Uq that comply with the condi¬ 
tion of Eq. ( |21a| ) and eventually make up the feasibility 
set have to lie somewhere in this conical volume. Now 
by inspection, for a given Uq, Eq. ( |21b ) states that the 
points in it lie in an inverted cone whose apex is ( 770 , vq). 
Thus, referring to Fig. [laj the question now is: does there 
exist a point ( 770 , Wo) in the light gray cone which is the 
apex of an inverted (dark gray) cone that includes all the 
points in it? 

If we consider Uq = 17(1,0) (heterodyne detection, 
or QSD), the two cones intersect each other right at 
\v\ = r] = 0.5; see Fig. lb It can obviously be seen 
that a set it containing only unravelings with efficiencies 
77 < 0.5 is not capable of demonstrating EPR-steering (as 
they are confined in the shared volume of the two cones). 
In other words, to be able to rule out all OPDMs, at 
least one efficiency has to be greater than 0.5. However, 
as mentioned in Ref. m this is just a necessary condition 
and is not sufficient to prove the detector dependency of 
stochastic conditional states g. This insufficiency can be 
seen immediately by looking at the upper half part of 
Fig. lb where there are other unravelings in the over¬ 
lapped volume that are also a coarse graining of 17(1,0), 
those with |t>| < 1 — 77 < 0.5. 



FIG. 2. (Color online). Steering parameter for K = 3 
measurement settings: the temporal behavior of the steering 
parameter S 111 , Eq. (231, for different values of the detection 
efficiency rj. It is evident that all unravelings for which 77 < 0.5 
are not able to demonstrate EPR-steerability. Nevertheless, 
once the system settles in the steady state (t —> 00 ), it still 
requires great improvements on detection efficiency to the ex¬ 
tent that it should roughly satisfy 77 > 0.8. Inset shows the 
value of S m versus 77 for the system in g ss . 


IV. DEMONSTRATING DETECTOR 
DEPENDENCE 

Having determined necessary conditions for being able 
to demonstrate detector-dependence using diffusive un¬ 
ravelings, we now turn to sufficient conditions. This will 
involve considering particular open qubit systems, cou¬ 
pled to one or more environments, with two or more dif¬ 
ferent unravelings, as well as considering various EPR- 
steering inequalities that might be violated. 


A. L = 3, I< = 3 

Consider first an open qubit system with L = 3 ir¬ 
reversible environmental channels, with the conditional 
evolution of system’s state characterized by the follow¬ 
ing SME 

3 

dg = y: (v[y/^ai\gdt + li[^aidVi\gy (22) 

l=i 


We consider the case where the complex Wiener pro¬ 
cesses are taken to be in the form of Eq. § with L = 3. 
Specifically, we consider performing three different un¬ 
ravelings, labelled by k £ {1,2,3}, in which the experi¬ 
menter obtains information about tJk from the by choos¬ 
ing (f>k = 0 for the £:th component of the observable, while 
giving rise to noisy Hamiltonian evolution for the other 
components by choosing (f>i = tt/2 for the l k. Every 
one of these three monitorings can be encapsulated in 
an unraveling matrix as shown in Eq. ([5]). For instance, 
for the k = 3 unraveling it is simplified to 0 = 77 /, and 
T = 77 diag(—1, — 1,1). The intuition behind this choice 
is that each unraveling yields information only about one 
observable, and so would be expected to make this ob¬ 
servable take a well-defined value. This would lead to 
quite different conditioned states for the different unrav¬ 
elings, which is what is needed to disprove all OPDMs. 

Thus to prove the detector-dependence of g a suit¬ 
able EPR-steering criterion should be examined. Such 
a criterion is often expressed in the form of an inequal¬ 
ity. For example, from the fact that the Bloch vector 
r = ((<Ji, < 72 , d 3 ) T ) must obey ||r|| 2 < 1, and the convex¬ 
ity of (&i) = Tr [gai\ in the conditional state matrix g , it 
can be shown that every OPDM satisfies mm 


Here < 7 ; are the usual Pauli operators. In comparison to 
the notation of Eq. ([ 2 ]), H = 0 in a relavent rotating 
frame, c = 02, <73), and dV = (dVi, dV-j, dV^. 

An unraveling equivalent to the one described by the 
SME given in Eq. (22) can be realized for a quantum 
system with one environment (for example, in a circuit 
QED setup) but by measuring the three Pauli observables 
concurrently mm- 


5 m = ^E fe [(d fe ) 2 ] <1. (23) 

fc=1 

Here E fc [*] means the ensemble average under the un¬ 
raveling k. That is, three unravelings are implemented 
to calculate the steering parameter 5 m . Note that due 
to the symmetrical nature of the problem the ensemble 
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FIG. 3. (Color online). Ensemble average of [(di ) 2 + (c^} 2 ]: 
results of numerical simulations of the first term in the L.H.S 
of the inequality given in Eq. (271 using X- and Y-homodyne 
detection as a function of time for r) = 0 . 6 . 


average does not depend on k when the system is in the 
steady state for t —> oo. 

Details of calculations are presented in App. |A l] Asa 
function of time t we plot in Fig. [2]the steering parameter 
given in Eq. ( [23] ) for various values of the detection effi¬ 
ciency 77 . It can be clearly seen that as long as 77 < 0.5 the 
EPR-steering inequality is never violated. This should 
not be very surprising as it was comprehensively ana¬ 
lyzed in the previous section. Nevetheless, violation oc¬ 
curs while the system is in the transient state provided 
that rj > 0.5. In contrast, once the system relaxes to its 
steady state (which might be much more interesting from 
practical points of view) only if r] > 0.8 it is possible to 
establish EPR-steerable states. This is also illustrated in 
the inset where the value of S ' 111 versus rj is drawn when 
the system dynamics reach to a stable point. The crit¬ 
ical efficiency of rj c « 0.82 is required to prove detector 
dependency of g in the steady state. 


B. L = 1, K m 2 

The high efficiency required by the 3-measurement set¬ 
ting suggests looking at the other extreme in which an 
open quantum system decoheres just to a single environ¬ 
ment. Consider that the SME of such a system is given 
by 


dg = 'D[ y /^a 3 ]gdt+ H[^ya 3 dV]g, (24) 

where again in notation of Eq. © Hamiltonian becomes 
zero in a suitable rotating frame, c = and the 

Wiener process is given by Eq. ^ with L = 1. Here 
we consider just two different measurement strategies, 
corresponding to (j) = 0 (X-homodyne, giving maximum 
information), and <j> = 7 r /2 (Y-homodyne, giving pure 


noise). Since L = 1, the correlation matrices are just 
scalars, 0 = rj = ±T for the two respective cases. 

For the 2-measurement settings we study a class of 
EPR-steering inequalities similar to the one given in 
Eq. (23) which in general can be formulated as 


E x 

where 


^ ^ ®-k{&k) 


k =1 


E 


Y 


E (1 - a k)(dkY 


k—1 


VJfe, 


0 < a k < 1 . 


<1, (25) 


(26) 


This shows that there are an infinite number of such in¬ 
equalities. However, one should carefully choose an op¬ 
timal one so that the left hand side of Eq. (25) consti¬ 
tutes a maximum. This is determined by whether the 
X-homodyne or Y-homodyne measurement scheme gives 
a larger value for the ensemble average of (dfc) 2 . If the 
former is the case one should choose a*, = 1 , and in 
the case of the latter ak = 0. There are protocols in 
which both X- and Y-homodyne detection give an equal 
amount of information about the fcth observable. In such 
a case ak £ [0,1] is arbitrary. Note that the comparison 
between E x [(<ffc) ] and E Y [(dfc) ] can be done through ei¬ 
ther numerical simulation or analytical calculation where 
possible. 

In this work we have chosen the optimal EPR-steering 
inequalities of the form (25) by analyzing the numerical 


results. It turns out that for the SME given in Eq. (24) 
the following inequality is optimal 

S 11 ee E y [{a,) 2 + <d 2 ) 2 ] + E x [(d 3 ) 2 ] < 1. (27) 

This can be clearly seen in Fig. [3] where the ensemble 
average of (di ) 2 + (d 2) 2 is plotted using both X- and Y- 
homodyne monitorings. The same argument holds for the 
second term of Eq. ( |27[ but this time unraveling using 
X-homodyne gives more information about <73 in com¬ 
parison to Y-homodyne. Each one of these two terms 
in Eq. (27) can be realized in an an EPR-steering ex¬ 


periment as follows. Bob randomly chooses between the 
two measurement settings, and so tells Alice to imple¬ 
ment either X- or Y- homodyne. Bob then correlates the 
outcomes of his tomographic measurements with those 
obtained by Alice. For details, see Ref. [3]. 

In Fig. [4] the time variation of the steering parameter 
S ' 11 for different values of detection efficiency 77 is shown. 
There are two key features closely resembling those of 
Sec. |IV A| Firstly, the necessary condition 77 > 0.5 for 
ruling out all OPDMS is satisfied. Secondly, the tran¬ 
sient phase of evolution establishes EPR-steering states 
provided the detection efficiency is greater than 0.5. 

However, in a third respect this case is quite different: 
in the long-time limit it is impossible to prove detector 
dependency of g no matter how close to unity 77 is. This is 
because the Y- homodyne unraveling generates just pure 
noise, which contains information about the system evo¬ 
lution (the noisy Hamiltonian) but no information about 
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FIG. 4. (Color online). Steering parameter for K — 2 
measurement settings: results of numerical simulations of the 
steering parameter S 11 , Eq. (271, as a function of time for 
different values of rj. Here for the sake of clarity we discard 
labels in between solid curves, but they follow the same color 
style and order as in Fig. [2] When the detection efficiencies 
satisfy r/ < 0.5 there is no possibilities (for homodyne unrav- 
elings) to demonstrate EPR-steering, as expected. However, 
beyond this limit, it is possible to observe EPR-steering if the 
system is in its transient phase of evolution. In contradiction 
to the case of L = 3, K = 3, in the long-time limit there is 
never any violation of Eq. (27|. Curves for r\ < 0.3 illustrate 
that the time taken to reach the asymptotic value of 1 grows 
as the efficiency gets smaller. 


the system state. Thus for any 77 < 1, the system purity 
decays deterministically over time under the <j> = 7t/ 2 
unraveling, eventually reaching a completely mixed state 
as t —> 00 . In this limit, the first term in Eq. (27) will 
be zero, and since the second term is bounded above by 
unity, there is no way the inequality can be violated. 
Violations occur at short times only because we started 
the system in a pure state. Further details are given in 
App. |A 2[ It is interesting to note that the transient evo¬ 
lution here spans over an interval several times longer 
than the duration in which S ' 111 attains steady state. 


C. L = 2, K = 2 

Finally, we concisely present the results of our studies 
of an open qubit system coupled to two irreversible chan¬ 
nels (L = 2), and compare them to the two results of us 
previously obtained in Ref. [15] . To this end consider a 
qubit SME in the form of 

dg = (v[y/vvi]e + K[VTmdVi]o), (28) 
z=± 

where a± = (dy ± i<T 2)/2 are raising and lowering opera¬ 
tors. As before, the system’s conditional state matrix g 
evolves in an appropriate rotating frame so that H = 0 


in the notation of Eq. ([ 2 ]). Also the Lindblad vector is 
c = (y/yZa-, yjy +<r+), and the independent Wiener pro¬ 
cesses are given by Eq. § Thus, depending on the 
measurement strategy, determined by the phase term (f> 
in Wiener process, the average current of the two differ¬ 
ent decoherence channels would be 

(J±)=VV 7l(<Vt+e i20 <7 T ), (29) 

such that the X-homodyne detection (<j> = 0) gives infor¬ 
mation about g 1 

(J±f=vV7 ±<*i>, (30) 

and the outcome of the Y-homodyne detection (</> = 7 t/2 ) 
is given by 


<J ± ) y = ±</7±(d 2 ). (31) 

This SME has been studied in Ref. m in order to pro¬ 
pose an experimental test for disproving all OPDMs us¬ 
ing quantum jumps unravelings for experimental setups 
with detection efficiencies less than 50%. In this protocol, 
however, Alice should be able to realize n measurement 
settings in the azimuthal direction <p € {(j/n.)7r}™ =1 and 
one in the z direction to be able to work out the consid¬ 
ered EPR-steering inequality m 


1 

n 


j =1 


1(^)1 


f{n)W 


\A ~ (^a ) 2 


< 0. (32) 


Here <j v = <j-e lLp + <J+e , and the function f(n) is 
given in Ref. m which monotonically decreases with n, 
asymptoting to 2/7T as n —> 00 . Thus in the limit of large 
n a smaller detection efficiency is required to demonstrate 
detector dependence of the conditional state of an open 
quantum system. It turns out that the critical efficiency 
is approximately r] c « 0.59 provided that one is capable 
of performing an infinite number of monitoring schemes, 
and in addition the condition R = y + /q_ -C 1 is satisfied. 
These assumptions are nontheless not very convenient. 
The former is not practical as it requires great number 
of measurement settings and the latter leads to a very 
small violation in comparison with the maximum possible 
violation when the detection is perfect and R = 1 m- 
However, it was shown that with just a finite number of 
measurement settings, K = n + 1 = 5, one would be 
able to have a decent (0.05) violation of Eq. (32) with 
?7 ss 0.78. In what follows we show that an even more 
suitable EPR-steering inequality can be derived to rule 
out all OPDMs with a sufficient condition as close as 
possible (from above) to the necessary condition bound 
of 50%. 

Thus for the case given by Eq. (28), because the a± are 
not Hermitian, both the X- and Y-homodyne measure¬ 
ments yield information about the system, but about dif¬ 
ferent observables, namely a\ and 02 respectively, given 
in Eqs. (30) and (31). This suggests yet another varia¬ 
tion derivable from the inequality ||r|| 2 < 1, namely the 





















following EPR-steering inequality 


si- s + <^] + E y <i 2 ) 2 +Ml 


] < 1, (33) 


which is a special case of Eq. (25) with a\ = 1 , 0:2 = 0 


and 0:3 = 1/2. This assignment should be now obvious 


considering the argument made in Sec. IV B and recalling 
Eqs. (30) and ( |3T] ). In fact the value of 013 is arbitrary 
as both X- and Y-homodyne give the same amount of 
information about < 73 . The choice a 3 = 1/2 emphasizes 
the symmetry. 

For sufficiently small R <C 1 it is easy to derive an 
analytical relation for the steering parameter in Eq. (331, 



S/t 1 =* 1 + 8(77 - 0.5)i?. 


(34) 


The details of the calculation are presented in App. [ 
The analytic relation, Eq. ( [34] ) , cries out not only the nec¬ 
essary condition mentioned in the Sec. |III[ but also shows 
that 77 > 0.5 is sufficient to rule out all OPDMs. How¬ 
ever, this result holds only for small R’s and as such it will 
not be convincible enough from experimental viewpoint 
due to a very small violation. Therefore, a desired ob¬ 
jective would be to achieve as low as possible to r] c ~ 0.5 
and obtain a decent violation of Eq. (33) so that realiz¬ 


ing experimentally EPR-steerable states, and hence To 
this end, for a decent violation of Eq. (331 an optimal 


value of R is obtained. For example, a few percent of 
the difference between the threshold and the maximum 
possible value allowed by quantum mechanics (here unity 
and 2 , respectively) should be almost enough to take into 
account imperfections. It turns out that an efficiency 
of approximately 0.72 is needed to acquire a 5% viola¬ 
tion with the ratio of transition rates being optimized 
at R ~ 0.2, see Fig. [Hj In this figure we plot S± as a 
function of 77 with the critical efficiency r] c = 0.68 for vi¬ 
olating Eq. (33). Also, the inset (bottom-right) depicts 
the same steering parameter versus R for a fixed detec¬ 
tion efficiency 77 = 0.72. 

Finally, one might question the validity of approxima¬ 
tions made to obtain Eq. (34), in the sense that it is 


implying that 77 > 0.5 should be sufficient for disproving 
all OPDMs in the limit where R is very small. To sup¬ 
port these simplifications which led us achieve the above 
mentioned analytic formula we present the outcome of 
numerical simulations of Eq. (33) for R ss 0.01 in the 
inset (top-left) of Fig. [HJ As it can be seen S// varies 
linearly with respect to 77 and this is exactly the behav¬ 
ior one would expect from Eq. (34). Moreover, this plot 
demonstrates that the critical efficiency is r) c ss 50.7% 
which is very close to what analytical calculation pre¬ 
dicts, that is 50%. This is a very good improvement on 
the critical efficiency required for proving detector de¬ 
pendence of stochastic conditional states g. 


FIG. 5. (Color online). Steering parameter for K = 2 
measurement settings and L = 2 de coh erence channels into 
the environment: we plot S±, Eq. (331, as a function of 77 
for an optimal value of R « 0.2 with a significant violation 
of 5% at rj = 0.72; see text for details. Critical efficiency 
here is r/ c = 0.68. Inset: top-left shows the same plot but for 
R « 0.01 with critical efficiency of r/c ~ 0.507; bottom-right 
displays variation of S± versus R for 77 = 0.72. 


V. CONCLUSION 

In summary, we have investigated the dependence of 
diffusive unravelings of open quantum systems on the 
existence of a distant detector. Due to imperfect detec¬ 
tion schemes some different measurement strategies on 
an open qubit system were studied to obtain a minimum 
quantum efficiency required to disprove objectiveness of 
pure-state evolutionary models. In this work we sought 
two objectives. Firstly, to determine some general nec¬ 
essary conditions for proving detector dependency of dy¬ 
namical quantum events.And secondly, to derive suffi¬ 
cient conditions by considering different protocols. 

We derive necessary conditions for being able to 
demonstrate EPR-steering in the form of a no-go for in¬ 
efficient diffusive unravelings. These conditions depend 
only on the parameters of the unravelings, not on any 
details of the system. This theorem can be applied in its 
most general form by recasting the condition as a special 
type of semidefinite programming known as a feasibility 
problem. The simple case of L = 1 can be thoroughly 
analyzed by a graphical representation involving right 
cones. 

In contrast to the case of necessary conditions, deriving 
sufficient conditions can only be done case by case, in that 
they depend on the systems being examined. To achieve 
this goal three factors may be inspected for an open qubit 
system: varieties of irreversible evolution equations, dif¬ 
ferent measurement schemes, and various EPR-steering 
inequalities. The threshold efficiencies required rj c for 
each scenario along with those of other previously pro¬ 
posed experimental tests (for diffusive unravelings) are 
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TABLE I. The critical efficiencies r/ c required for violating 
different EPR-steering criteria. Except for the last two rows, 
the rest represent entirely different physical systems coupled 
to different number of decoherence channels; dissipative (Diss) 
and non-dissipative (ND). Also for each case specific measure¬ 
ment schemes are considered. 


Decoherence Measurement Steering 

channels (L ) settings (K ) parameter 


3 (ND) 3 

1 (ND) 2 

1 (Diss) 2 

2 (Diss) 2 

2 (Diss) 2 


S m , Eq. <[23]> 

S’", Eq. pf 
Eq.(3) of Rirg| 
S%, Eq. {33} 
Eq.(8) of Ret. 1151 


->■ 00 
Vo _ 

0.82 

0.73 

0.5 

0.59 


summarized in Table [TJ The investigated environmental 
decoherence channels fall into two categories; dissipateve 
(Diss) and non-dissipative (ND). According to the re¬ 
sults obtained in this work, the former class shows better 
figures for the critical efficiency. Among all of these the 
case in which the system decoheres into two environments 
and the latter undergo monitoring via two measurement 
settings (Sec. IV C) shows a promising candidate for be¬ 
ing probed empirically. Lastly, this work may, on the 
one hand, encourage experimentalists to endeavour to 
increase the efficiencies of detection schemes, and, on the 
other hand, stimulate theorists to invent still more robust 
EPR-steering tests. 
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Appendix A: Simulating averages 

1. L = 3, K = 3 

For the three measurement settings general dyne un¬ 
raveling as mentioned in the main text, depending on the 
measurement strategy, Eq. ([ 5 ]) is written in the following 
forms 


U \ k=1 = 77 diag(l, 0,0,0,1,1), (A.l) 

U \ k=2 = 77 diag(0,1,0,1,0,1), (A.2) 

U\ k= 3 = 77 diag( 0 , 0 , 1 , 1 , 1 , 0 ), (A.3) 

which are clearly PSD. Then the conditional state of the 
system which evolves according to Eq. ( [22] ) can be rewrit¬ 
ten in the Bloch representation so that dynamics of the 
system is governed by the following nrulti-dimensional 
stochastic differential equation (MDSDE) 

dr = A(r)dt + B(r)dW (f), (A.4) 


where 


r = (x, y, z) T = ((cti) , (<f 2 ), (d 3 )) T , (A.5) 

W = (1T 1 ,1T 2 ,1T3) T , (A. 6 ) 

A = —4r, (A.7) 

and the WjS are independent Wiener processes. The ma¬ 
trix B also takes the following form 


B = 2 y/rj 


0 —2 —xz 

z 0 — yz 
—y x 1 — z 2 


(A. 8 ) 


Here, we have just expressed details of Eq. (A.4) for the 


k = 3 case. Note that the ultimate goal is to calculate 
the steering parameter S U1 , given in Eq. (23), composed 


of ensemble averages. By virtue of symmetrical nature 
of the problem the ensemble average E fc [(ttfe) 2 ] is inde¬ 
pendent of k when the system evolves, for a sufficiently 
long time. Also note that the fcth measurement strat¬ 
egy merely gives noisy output for the rest of observables. 
That is, E fc [(iTfc') 2 ] = 0 for k' 7 ^ k. 

As it was discussed in Sec. |IV A the kth unravelling of 
the SME, Eq. (22), gives only information about the 07 


observable. Therefore, to calculate E fc [(o 7 ) 2 ] in Eq. (23) 
it is convenient to introduce some variables /3 ,6 and 7 
with the following relations 

(<j m ) +i{a n ) = \fpe l6 , for m,n^k (A.9) 


( d - fc ) 2 = 7. 


(A.10) 


This allows us to work with a MDSDE which is more 


relevant to what we need to work out Eq. (23). The 


following lists different matrices that characterize such a 
MDSDE for k = 3 


r=(/3,0, 7 ) T , 

W = (wp, Wg, H / 7 ) T , 


and 


A = 


8777 + 4 t 7/?7 — 8/3 
^ 777 2 sin(40) 

477/3 + 477(1 — 7) 2 — 87 


(A.ll) 

(A.12) 


(A.13) 


B = 4^77 


-VWi 0 —0y/l 

0 vr 7m 0 

VWy 0 ^7(1 - 7) 


.(A.14) 


Therefore, to evaluate S ' 111 in Eq. ( |23| the ensemble 
average of 7 for each individual unraveling, E fe [(dfc) 2 ], 
should be calculated. As it follows from Eqs. (A.ll) 


(A.14), the evolution of 7 in turn is coupled to that of 


/3 and is independent of variations of 9. Thus, a set of 
just two coupled SDEs for j3 and 7 is simulated. This 
MDSDE is solved numerically using the Milstein method 
[32] . For sufficiently long time that the system is in the 
steady state data are sampled and the ensemble averages 
are calculated. 
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2 . L = 1, K = 2 

The unraveling matrix based upon homodyne measure¬ 
ments on the output held of the system described by 
Eq. (24) becomes U = ?ydiag(l,0). Then th e sto c hastic 
equations of motion in the notation of Eqs. (A.4)- (A.6) 
can be expressed as 


(A.15) 



' 1' 


0 0 — xz 

A = -2 

1 

, B = 2^/77 

0 0 —yz 


0 


0 0 1-z 2 


Now using Eqs. (A.9)- (A.12), this MDSDE is trans¬ 
formed into 


A = 4r; 


B = A^rp/ 


7/3 - f3/r] 

|sin(4 9) 

(1-7) 2 

0 0 

0 l/y/Af3 0 

0 0 1-7 


(A.16) 

(A.17) 


The same argument as the three measurement settings 
holds here where the system dynamics at any given time 
can be described by variables (3(t ) and 7 (f), and also the 
same procedure is employed to calculate the ensemble 
average of [(<t 3 ) 2 ] ■ 

For the noisy Hamiltonian evolution measurement 
scheme though, the unraveling matrix is U = 77 diag(0,1). 
The ensemble average of [( ay ) 2 + ( 02 ) ], which is equal 
to that of [/?] through Eq. (A.9), is obtained by solving 
the following SDEs 


d/3 = 4/3 (?/ — 1) dt, 
d'y = 0 , 

which lead to the equations of motion 

7 (t) =7(0)- 


3. L = 2, K = 2 


(A.18) 

(A.19) 


(A.20) 
(A.21) 


following correlation matrices 

0 = diag(? 7 , 77 ), T = ?ydiag(e -21 ^, e 2l< ^). (A.22) 

These diffusive unravelings yield information only about 
observables <5y (/ i> = 0) and 02 (</> = 7 t/ 2 ). In fact, one of 
the reasons why we are working with the specific EPR- 
steering inequality given in Eq. (33) is because there is 


not an appropriate unraveling to be useful for calculating 
the ensemble average of (a 3 ) . However, as it is the case 
here, one could still estimate it from the measurement 
record for ay and & 2 - 

Due to the symmetry of the problem both terms in 
the L.H.S of Eq. (33) return the same value. This means 
we can evaluate one and double it. Thus, the case of X- 
homodyne with (f> = 0 is studied. The dynamics of the 
qubit is confined to x-z plane such that at any instant of 
time its state can be described by a point ( x(t ), 0 , z(t )) in 
the Bloch sphere, where the Bloch vector is r = ( x , y, z). 
The trajectory of this point in the sphere is tracked by 
solving the coupled SDEs 

dx = —(te/2) xdt + y/rpyZ( 1 + z — x 2 )d,W- + 

v^7+(l - 2 - x 2 )dW+, (A.23) 

dz = (— 7 s ^ + 7 a )dt - y/rpfZx{l + z)dW- + 

\/Wi+x{ 1 - z)dW+, (A.24) 

where 7 e = 7+ + 7- and 7 a = 7+ — 7^ . As we have 
done for previous simulating jobs we use Milstein method 
to numerically compute these equations. Once the sys¬ 
tem reaches to its steady state, t —> 00 , data is being 
recored to calculate E x [(iT 2) 2 + 1 / 2 (<t 3 ) 2 ], and hence 
E Y [(dy ) 2 + 1/2(<t 3 ) 2 ], for each measurement scheme as 
described above. 

Finally, when R = 7+/7_ <C 1, the conditional state 
of the system under diffusive unraveling is almost in 
proximity of the ground state so that approxim ations 
x = O(Vr ) and 1 + z = O(R) hold. Then Eqs. (A.23) 
and ( |A.24[ ) are transformed, to leading order, to 

dx = — ( 7 / 2 ) x dt + 2 \JrjjR dW +, (A.25) 

dz = "f[2R — z — l]dt + 2^/rp/R.x dW + , (A.26) 

where 7 = 7-. It is straightforward to calculate the 
ensemble averages E[x 2 ] ss AyR and E [z 2 ] ~ 1—AR for the 


steady state. Plugging these results back into Eq. (33) it 


In this scenario the homodyne detections (with phases 
and —(j>) of the two environmental channels gives the 


is easy to show that for small R it is approximated by 
S R<< 1 “ 1 + 8 ( 7 / - 0.5)33. (A.27) 
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